Quantum many-body theory of qubit decoherence in a finite-size spin bath 
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Decoherence of a center spin or qubit in a spin bath is essentially determined by the many-body bath evolu- 
tion. We develop a cluster-correlation expansion (CCE) theory for the spin bath dynamics relevant to the qubit 
decoherence problem. A cluster correlation term is recursively denned as the evolution of a group of bath spins 
divided by the cluster correlations of all the subgroups. The so-defined correlation accounts for the authentic 
(non-factorizable) collective excitations within a given group. The bath propagator is the product of all possible 
cluster correlation terms. For a finite-time evolution as in the qubit decoherence problem, a convergent result 
can be obtained by truncating the expansion up to a certain cluster size. The two-spin cluster truncation of 
the CCE corresponds to the pair-correlation approximation developed previously [Phys. Rev. B 74, 195301 
(2006)]. In terms of the standard linked cluster expansion, a cluster correlation term is the infinite summation of 
all the connected diagrams with all and only the spins in the group flip-flopped, and thus the expansion is exact 
whenever converges. When the individual contribution of each higher-order correlation term to the decoher- 
ence is small (while all the terms combined in product could still contribute substantially), as the usual case for 
relatively large baths where the decoherence could complete well within the bath spin flip-flop time, the CCE 
coincides with the cluster expansion [Phys. Rev. B 74, 035322 (2006)]. For small baths, however, the qubit 
decoherence may not complete within the bath spin flip-flop timescale and thus individual higher-order cluster 
correlations could grow significant. In such cases, only the CCE converges to the exact coherent dynamics of 
multi-spin clusters. We check the accuracy of the CCE in an exactly solvable spin-chain model. 
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I. INTRODUCTION 

The decoherence of a center spin in a spin bath 1 has 
been of interest in spin resonance spectroscopy for a long 
histor y 2 ' 3 ' 4 ' 5 and is also a paradigmatic model in studying 
the state collapse in quantum mechanics. Recent revisited 
interest in this problem is mostly due to the decoherence 
issue in quantum computing. 6,7,8 Most relevant are single 
electron spins in quantum dots or impurity centers where 
the dominating decoherence mechanism at low temperatures 
(such as below a few Kelvins) is the nuclear spins of the host 

1 tt j 9.10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34 

In a small system such as a quantum dot, the center spin 
(hereafter referred to as qubit for clarity) and the spin bath, in 
the time-scale of decoherence, is a relatively isolated subsys- 
tem in the whole environment. Thus the qubit decoherence 
is due to the entanglement with the bath during the coherent 
evolution of the whole system^ 

In this paper, we are interested in the so-called pure dephas- 
ing in which the qubit experiences no longitudinal relaxation 
but only loses its off-diagonal phase coherence. The pure de- 
phasing is relevant to an electron spin under a moderate or 
strong magnetic field (> 0.1 T for a typical GaAs dot), where 
the electron spin flip due to the hyperfine interaction is largely 
suppressed by the Zeeman energy mismatch between the elec- 
tron and the nuclei. In the absence of qubit flip, a qubit-bath 
system has a Hamiltonian as 

H = \+) #«<+| + |-> H^(-\, (1) 

by which the spin bath is driven by different Hamiltonians 
i/ (±) depending on the qubit states |+). When a coherent qubit 
state C+ |+) + C |-) is prepared, the initial state of the qubit- 
bath system is the product state (C + |+}+C_ At time 
T, the bath evolution \J) U^T)) = e - iH{± ' T \J) predicated 



on the qubit state |±) establishes an entangled state C + |+) ® 
\J + (T)) + C |-) ® \J~(T)). The qubit coherence is reduced 
from P+ _(0) = C + C! to p + _(T) = C + C! (J-(T)\J + (T)). 
The decoherence is characterized by the bath state overlap 
(^T (T)\^T + (T)). Thus the key is the many-body bath dynam- 
ics caused by the interaction within the bath. Without many- 
body interactions, the bath would not evolve (except for a 
trivial phase factor that can be eliminated by standard spin- 
echc 35 ) and the qubit coherence would not decay. 

Recently, a variety of quantum many-body theories for 
nuclear spin bath dynamics have been developed including 
the density matrix cluster expansion (CE) , 17 ' 18 ' 19 the pair- 
correlation approximation , 20 ' 2 1 ' 22 and the linked-cluster ex- 
pansion (LCE). 23 In the pair-correlation approximation, each 
pair-wise flip-flop of nuclear spins is identified as an elemen- 
tary excitation mode and is taken as independent of each other. 
To study the higher oder correlations, the Feynman diagram 
LCE is developed. The evaluation of higher-order LCE, how- 
ever, is rather tedious due to the increasing number of dia- 
grams, especially for spins higher than 1/2 (see the Appendix 
for details). The density matrix CE is developed in the spirit 
of the standard cluster expansion or virial expansion for in- 
teracting gases in grand canonical ensembles. 36 ' 37 ' 38 It serves 
as a simple method (without the need to count or evaluate 
Feynman diagrams) to include the higher-order spin interac- 
tion effects beyond the pair-correlation approximation. The 
CE calculations show that when the pair-correlation contri- 
bution to the decoherence is suppressed by pulse control, 21 
the higher-order correlations are not negligible. 19 In order to 
obtain the decoherence exponent from the CE, however, the 
terms involving overlapping clusters have to be neglected, 
which make the accuracy problematic even when the expan- 
sion converges. This problem limits the CE to applications in 
large baths in which the contribution of each individual clus- 
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ter to the decoherence is small and the overlapping correction 
is unimportant. 

In this work, we develop a cluster-correlation expansion 
(CCE) method in which the bath spin evolution are factor- 
ized into cluster correlations. Each cluster correlation term 
is equivalent to the sum of all the LCE series^ 3 - consisting 
of a given set of bath spins flip-flopped. In particular, the 
two-spin cluster correlations include all diagrams with two 
spins flip-flopped and is equivalent to the pair-correlation 
approximation. 20,21,22 The CCE bears the accuracy of the LCE 
(the results are accurate whenever converge) and the simplic- 
ity of the CE (without the need to count or evaluate Feyn- 
man diagrams), while free from the large-bath restriction of 
the CE. The CCE coincides with the CE in the leading order 
of the short-time expansion, which is applicable in large spin 
baths where the decoherence completes well within the bath 
spin flip-flop time. For small baths, however, the qubit deco- 
herence may not complete within the bath spin flip-flop time 
and individual higher-order cluster correlations could grow 
significant. In this case only the CCE converges to the exact 
coherent dynamics of multi-spin clusters. Such coherent dy- 
namics of small clusters of bath spins is of special interest in 
systems with randomized qubit-bath couplings. An interest- 
ing example is nitrogen-vacancy centers in diamond s 39 ' 40 ' 41 ' 42 
which are coupled to randomly located bath spins (carbon- 13 
and nitrogen nuclear spins) in the proximity. 

This paper is organized as follows. In Section II, we derive, 
for a generic spin bath Hamiltonian, the CCE from a recur- 
sive cluster factorization procedure. We show that the CCE 
is equivalent to an infinite resummation of the LCE series and 
also compare it to the CE. In section III, we check the accuracy 
and convergence of the truncated CCE in an exactly solvable 
model (the one-dimensional spin-1/2 XY model). Section IV 
draws the conclusions. The Appendix contains the details of 
the LCE for spin bath dynamics, extended to spins higher than 
1/2. 



II. CLUSTER-CORRELATION EXPANSION 
A. Motivation: Pair-correlation approximation and beyond 

As discussed in the Introduction, the Hamiltonian for the 
pure dephasing problem has the form of Eq. £[). For a given 
initial bath state U) (which could be considered as one sam- 
ple chosen from a thermal ensemble), the qubit coherence is 
characterized by 



m->T -«/<+> r 



£.(?) = <J\e M 'e 



\J). 



For a thermal ensemble of baths, a further ensemble average 
should be processed. The thermal fluctuation leads to the in- 
homogeneous broadening, which can be eliminated by spin 
echo. To focus on the qubit decoherence due to the quantum 
dynamics of the bath, throughout this paper, we consider the 
decoherence for a single bath state without the ensemble aver- 
age. For temperatures much higher than the bath spin flip-flop 
rates (~1(T 9 K for nuclear spins in GaAs), the thermal ensem- 
ble has no off-diagonal coherence and U) can be taken as a 
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FIG. 1 : Visualization of a cluster containing (a) two or (b) three bath 
spins (black arrows). The spins outside the cluster (gray arrows) are 
taken as frozen in calculating the cluster contribution. 



noninteracting product state U) = ®„ \j„), where j„ denotes 
the quantum number of the nth bath spin quantized along the 
external magnetic field. 

In this subsection we illustrate the central idea of the CCE 
method (as an extension of the previously developed pair- 
correlation approximatio n 20 ' 2 1 ' 22 ) using a bath consisting of 
N spins Ji, J2, • • • , and Jjy with only pairwise secular interac- 
tions 

ff (±) =±l+E( z » ± i) / « + zK» ± % i )« 

n nvtn * ' 

+ ^ l#m,n ± ) J m J n> 



(2) 



where Q. (Z„) is the qubit (bath spin) splitting energy, z„ is 
the diagonal qubit-bath spin interaction constant (correspond- 
ing to the hyperfine interaction strength in an electron-nuclear 
spin system), Z\ m (B,^ ln ) is the diagonal (off-diagonal) in- 
trinsic bath interaction strength, and d n , m (b n , m ) is the diago- 
nal (off-diagonal) extrinsic bath interaction depending on the 
qubit states (which could result from the interaction mediated 
by virtual flips of the qubit spin while real flips are suppressed 
by the large energy mismatch 20 ). 

In the pair-correlation approximation, the qubit coherence 
is given by the product of all possible spin pair contributions 
up to a phase factor^ . 



n 

('■../■) 



(3) 



where the contribution due to the flip-flops of a spin pair [i, j] 
[see Fig. Q2 a)] is 



HJ) 



rr\ ih { r \T -ih)*\T\ cr\ 

J\e M e m \ J), 



(4) 



with the Hamiltonian governing the pair dynamics obtained 
from the full Hamiltonian i/ (±) by excluding the flip-flops of 
all spins other than J, and J ; 
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which is equivalent to replacing in the full Hamiltonian the 
spins outside the pair with their mean-field averages, 



Jj,{{j\3*u\ J)}). 



(6) 



The pair-correlation approximation is valid for situations 
where pair correlations dominate. When the collective flip- 
flops of more spins become important, we need to considei 
the higher-order correlation correction Xcon defined by 



V (U) 



(7) 



To illustrate how the high-order correction could be evaluated, 
let us consider a spin bath of only three spins {1,2, 3). Obvi- 
ously, the result is 



-£corr — 



L 



-U,2)^{2,3|M3,l) 



(8) 



This result motivates a definition of non-factorizable three- 
spin correlations due to collective flip-flops. For a three-spin 
cluster [i, j, k] in a bath [see Fig.QJb)], when all the spins out- 
side the cluster are frozen, the qubit coherence is 



I rr\ thP-u? -ih\ + \j\ cr\ 



with the cluster Hamiltonian 



+ ^ l B m,n ± ~7p J Jm J n 

m,ne[i,j,k\ 

obtained from the full Hamiltonian // (±) by replacing the spins 
outside the cluster with their mean-field averages. The au- 
thentic (or non-factorizable) three-spin correlation is singled 
out by excluding all the pair-correlations 
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If all such three-spin correlations are picked up, the qubit co- 
herence is given by (up to a global phase factor) 

£ ~ f| £{«,./) • Y\ hum- 

W,j] VJM 

Thus a systematic cluster correlation expansion is motivated. 



B. General formalism of cluster-correlation expansion 

We consider a generic bath Hamiltonian for the pure de- 
phasing problem 

#(±) =H (±) (h,h,--- ,Jw), 
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FIG. 2: (a) Visualization of an arbitrary sequence of controlling n 
pulses at tifo,' • • , (indicated by vertical lines) and the corresponding 
contour-time-dependent Hamiltonian for electron spin decoherence. 
(b), (c), and (d) exemplify the cases of free-induction decay, Hahn 
echo, and Carr-Purcell echo, respectively. 



which need not contain only pairwise interactions, or con- 
serve the spin angular momentum along any direction. For 
instance, multi-spin interaction terms like JjJjJ^ and non- 
secular terms like JfJ 1 - could be present. 

For qubit decoherence under the control of an arbitrary se- 
quence of 7r-pulses applied at t\, t%, ■ ■ • , as shown schemati- 
cally in Fig. HJa), the bath evolution predicated on the qubit 
state is given by 



j*(D> = u {±) \j) , 



where 



f/(±) = . . . e -'H l±] (ti-t 2 ) e -iH m (t 2 - ti ) e -iH {±) ti 

The qubit coherence at the end of the evolution is 



£ = (J 



J 



It can be written in the contour time-ordered form as 23 



(9) 



(10) 



(11) 



where ,% is the time-ordering operator along the contour 
C : — > T — > 0. As illustrated in Fig. Ufa), the contour 
time-dependent Hamiltonian H(f) alternatively switches be- 
tween // (+) and each time the qubit state is flipped by a 
7r-pulse or when the time direction is reverted at T. The ex- 
amples for free-induction decay, single-pulse Hahn echo, and 
Carr-Purcell echo are illustrated in Figs. [2] (b), (c), and (d), 
respectively. 

Following the idea illustrated in the previous subsection, 
the cluster correlations are recursively defined as follows. 

1 . The empty-cluster correlation 
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is a pure phase factor obtained from Eq. (fTTT > by replac- 
ing the bath Hamiltonian H(t) with its mean-field aver- 
age(J\H(t)\J). 

2. The single-spin correlation 

L{i) = L{i)/Le,, 

where 

L {i} = (j\$ c e- i i h ^ t)d <\j\ 

is obtained from Eq. (fTTT i by replacing the bath Hamil- 
tonian H({J n }, t) with 

%M = #(Ji,{cm*ij)},*). 

in which all the spin operators except J, are mean-field 
averaged. 

3. The two-spin (pair) correlation 

L{',Ji = L V,hl (^-zL { i}L{j)) , 

where 

is obtained from Eq. (fTTT i by replacing the bath Hamil- 
tonian H({J„), t) with 

AM© = H (h, Jj, {{J \J,mu}\j}} , t) , 

in which all the spin operators except J,- and Jj are 
mean-field averaged. 

4. So on and so forth, the cluster correlation for an arbi- 
trary set of bath spins C is defined as 

Lc = , (12) 

n V 

C'cC 

where 

*c = (j'|#*"^* cW *|j') (13) 

is obtained from Eq. (fTTT i by replacing the bath Hamil- 
tonian H({J n }, t) with 

he® = H ({J„ ec ) , {{j \intc\ J)},t), (14) 

in which all the spin operators outside the cluster are 
mean-field averaged or their flip-flops are frozen. 

Thus, by definition, the qubit coherence is factorized into 
all possible cluster correlations as 

x= n Lc - (i5) 

Cc{l,2,-,iV} 



Calculating the CCE to the maximum order L\\2.--.n) amounts 
to solving the exact bath propagator, which is in general not 
possible. In the decoherence problem, we consider a finite- 
time evolution and it often suffices to truncate the expansion 
by keeping cluster correlations up to a certain size M, as the 
Mth-order truncation of the CCE (M-CCE for short), 

£ iM) = \\ Lc, (16) 

\C\<M 

where \C\ is the number of spins contained in the cluster C. 
As an example, for the pairwise Hamiltonian in Eq. the 
lowest nontrivial order of truncation is 

£ (2) = U \\ = U \\ (L {iJ) /L ) , (17) 
{;,;'} 

which is the pair-correlation approximation. 

C. Relation to linked-cluster expansion 

Saikin et al. have recently developed an LCE method for 
the qubit decoherence in a spin-1/2 bath. 23 The detailed de- 
scriptions of the LCE for a generic spin bath are given in the 
Appendix. In general, the bath evolution can be factorized 
using Feynman diagrams so that 

£ = exp(7r), (18) 

where 

n=(j\&ce- i t° md, \j) , (19) 

\ ' ' / connected 

is the sum of all connected Feynman diagrams obtained by us- 
ing Wick's theorem on the series expansion of Eq. (fTTT >. Some 
of the connected diagrams up to the 4th-order have been eval- 
uated in Ref. ll23ll for free-induction decay and single-pulse 
Hahn echo with a spin-1/2 bath Hamiltonian in the form of 
Eq. The complexity for diagram counting and evalua- 
tion increases dramatically when considering higher-order di- 
agrams (see Fig. QT] in Appendix IA 2t or higher spins [see 
Fig.QjJb) in Appendix IA 3l for the second-order diagrams for 
a spin-1 bath]. 

We notice that each diagram can be expanded as the sum of 
diagrams involving the flip-flops of different clusters of spins. 
As an example shown in Fig. EJa), a third-order diagram in- 
volving the flip-flop of a spin pair contains diagrams for spin 
clusters (1, 2), (1, 3), • • • , where the numbers stand for the in- 
dices of the spins flip-flopped (i.e., the J z quantum number 
changed). Thus all the connected diagrams can be classified 
according to the spin clusters instead of the interaction orders. 
For an arbitrary cluster C, we define n(C) as the sum of all 
connected diagrams in which all and only the spins in clus- 
ter C have been flip-flopped. For instance, some of the dia- 
grams constituting 7r(0) and j) for a spin-1/2 Hamiltonian 
in Eq. (f2]i are shown in Figs. [3lb) and (c), respectively. With 
these {n} functions, the LCE is expressed as 

cc{i,2,-,jv) 
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(b) 7T(0) = 



(c) 7i{i,j) - 



+ v 




FIG. 3: (a) Expansion of a third-order connected diagram into dia- 
grams involving the flip-flops of different clusters of spins, (b) and 
(c) show the diagrams contained in jf(0) and n(i, j), respectively. In 
the diagrams, a solid arrow denotes the propagation of a spin. A 
wavy (dotted) line connecting two solid arrows/spins denotes a pair- 
wise off-diagonal (diagonal) interaction, and an open-ended line rep- 
resents the interaction with the qubit spin or an external field. Each 
filled circle/empty circle/empty square on a solid arrow denotes a 
J + /J~/J z operator, which raises/lowers/keeps invariant the J z quan- 
tum number of that spin. 



In particular, the infinite summation of all the connected dia- 
grams for a certain cluster C and all its subsets 



7T(C)= YjTtiC), 



(21) 



C'CC 



can be obtained from the series expansion Eq. ( fT9b by drop- 
ping all the terms involving the flip-flop of spins outside the 
cluster C, or, equivalently, by reducing the bath Hamiltonian 
H{t) to the cluster Hamiltonian /ic(0 in which the spins out- 
side the cluster are mean-field averaged. Thus we have 



e niO = j"! /(C) = Ijfa-'L j\ = U 



C'QC 



Comparing this to Eqs. <fT~2t - (TT~4-b , we immediately have 



L c = e 



MP) 



(22) 



(23) 



i.e., a cluster correlation term corresponds to the infinite par- 
tial summation of all the connected diagrams in which all and 
only the spins in the cluster have been flip-flopped. 

From the LCE expression of the CCE in Eq. ( l23b . it is obvi- 
ous that the short-time profile of the decoherence due to clus- 
ters of a certain size is determined by the lowest-order dia- 
gram. In particular, for free-induction decay with the secular 
pair-interaction Hamiltonian in Eq. (|2j, the leading order con- 
tribution from a cluster of size M is 



7f(C) ~ B M T M , 



(24) 



where B is the typical magnitude of the pair flip-flop interac- 
tion strength B m n + b mfl . If each spin interacts, on average, 
with q spins, then the number of size-M clusters is ~ Nq M ~ l , 
with N the total number of bath spins. The sum of all the 
leading order M-spin connected diagrams is ~ q~ l N ■ (qBT) M . 



Thus for qBT <k 1, i.e., for a time T much shorter than the 
bath spin flip-flop timescale 1/B, a truncated CCE converges. 
The short time condition T <K B l is usually satisfied for elec- 
tron spin decoherence caused by nuclear spins in typical quan- 
tum dots. The convergence of the truncated CCE, however, 
could go well beyond the short time restriction. One such sce- 
nario is small spin baths with disorder in qubit-bath couplings 
or in spin splitting energies, in which multi-spin correlation 
could develop at a time well beyond the short-time limit but 
the size of the contributing clusters could remain bounded due 
to the localization effect in a disordered system, as will be ver- 
ified later in this paper by numerical simulations. 



D. Relation to cluster expansion 

Witzel et al. recently developed a density matrix CE ap- 
proach to solving the nuclear spin dynamics in the electron 
spin decoherence problemists in the spirit of the cluster 
or virial expansion for interacting gases in grand canonical 
ensembles. 36 ' 37 - 38 Below we reproduce the basic procedure of 
the CE and compare it to the CCE. Instead of the ensemble 
CEinRef. [17], we consider a single sample state of the bath 
for a direct comparison. Defining the qubit decoherence due 
to a cluster C of bath spins as 



W(C) = 



\Lc\: 



(25) 



a hierarchy of cluster terms { W(C)} are recursively defined as 

W(i) =W(i), (26a) 
W(i, j) =W(i, j) + W(i)W(j), (26b) 
W(i, j, k) =W(i, j, k) + W(i)W(j)W(k) + W(i, j)W(k) 

+ W(i, k)W(j) + W(i)W(j, k), (26c) 



w(o =w(o + 2 n W(c ^ 



(26d) 



iCj) Cj 



where in the last line the sum runs over all possible parti- 
tions of the cluster C into nonoverlapping nonempty subsets 
Ci , C 2 , • • • . The Mth-order truncated CE (M-CE for short) is 



(27) 



{Cj}pj\<M C; 



with the sum running over all possible partitions of the bath 
into nonoverlapping nonempty clusters C\,C2, • ■ ■ of size up 
to M. 

In the cluster expansion for interacting gases in grand 
canonical ensembles with translational symmetry, the contri- 
bution from different clusters can be factorize d 36 ' 37 ' 38 and the 
evaluation of a truncated CE amounts to the calculation of a fi- 
nite number of finite-size cluster contributions, similar to the 
CCE in this paper. For a finite-size spin bath or for a bath 
without translational symmetry, however, such factorization 
of different clusters does not exist, which makes it essentially 
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impossible to calculate the sum in Eq. (l27l i even for a small- 
M-CE. For example, for a bath of N spins, the number of 
all terms containing only pair clusters is O(Nll) and all such 
terms have to be individually calculated and summed in the 
2-CE, which is practically impossible. 

A remedy is possible when all the cluster terms W(C) are in- 
dividually small (but the sum could still be substantial). Under 
such a condition, the CE can be approximated by a factorized 
form by adding some overlapping terms which are higher- 
order small. 17 For example, with the secular pair-interaction 
Hamiltonian in Eq. ©, for which W(i) = W(i) = 1, the M-CE 
is 

W (M) * ]~] [ 1 + . (28) 

1<\C\<M 

under the small-term condition 

\W(C)\ « 1, for |C| > 1. (29) 

Comparing the factorization approximation in Eq. ( l28l i to the 
exact M-CE in Eq. ( 1271 ). the error added is 

6W (M> = £ k)+ j, k)W(k, I) + ■ ■ ■ , 

/'< j<k i< j<k<l 

(30) 

containing products of any set of cluster terms sharing at 
least one spin, i.e., the overlapping terms. Such overlapping 
terms are higher-order small if each individual cluster term 
(for |C| > 1) is small. Furthermore, under the small-term con- 
dition, 1 + W(C) ~ exp [W(C)], and the CE for an arbitrary 
cluster becomes 

W(C) * Y] e^ c '\ (31) 

C'QC 

which takes the same form as the CCE. 

Thus the CCE coincides with the CE under the small-term 
condition Eq. ( |29l , which is justified for large spin baths 
where the number of contributing clusters is large and hence 
the contribution from each individual cluster remains small 
within the timescale of decoherence. The problem with the 
neglected overlapping terms is relevant for small spin baths 
where the coherent dynamics of a small number of multi-spin 
clusters dominating the decoherence may persist well beyond 
the bath spin flip-flop time and the small-term condition is no 
longer satisfied. In this case the CE will not converge to the 
exact multi-spin cluster dynamics, as will be seen in the nu- 
merical check in the next section. 



III. NUMERICAL CHECK 

Here we consider an exactly solvable spin bath model (the 
one-dimensional spin- 1/2 XY model) and compare the exact 
solution to the results obtained with the CCE and the CE. The 
TV-spin bath Hamiltonian conditioned on the qubit state |+) 
(with spin splitting constants dropped) is 

N N-l , , v 

fl« = ± X jJ z „ + X k ± y {j; +l r„ + J + n r n+1 ) , (32) 

«=l z «=l \ z ' 




Pulse delay time % 



FIG. 4: (Color online) The magnitude of size-ra cluster contributions 
(denoted as CCEm) to the exponential decoherence factor - In |X| in 
(a) free-induction decay (FID), (b) Hahn echo, and (c) Carr-Purcell 
echo (CP) for a long "sinusoidal" spin chain with N = 1000 spins. A 
solid (dashed) line indicates that the cluster contribution to - In |£| is 
positive (negative). The exact - In |_£| is also shown (empty squares) 
for comparison. 



where z„ denotes the qubit-bath spin interaction coefficient, B„ 
is the intrinsic bath interaction strength, and b„ is the interac- 
tion dependent on the qubit state. The initial bath state \ST) is 
taken as a product state of all bath spins, in which the orienta- 
tion of each bath spin is randomly chosen as up or down. The 
qubit-bath interaction coefficients \z n ) are either taken from 
a sinusoidal distribution z n = Zmax sin.(rm/N) (referred to as 
a "sinusoidal" spin chain) or randomly chosen from [0, z max ] 
(referred to as a "random" spin chain). Hereafter z max is taken 
as the unit of energy. The spin-flip interaction strengths \B n } 
and \b n ) are randomly chosen from [1(T 3 ,2 x 1CF 3 ], corre- 
sponding to typical bath spin flip-flop time r sl ~ 10 3 . The ex- 
act solution is obtained by the Jordan- Wigner transformation, 
which transforms the interacting spin- 1/2 chain to a noninter- 
acting fermion system42di 



A. Large spin bath 

As the first example, we consider a sinusoisal chain of 
TV = 1000 spins. In this case, the energy cost \z„ - z„+il f° r a 
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Pulse delay time x 

FIG. 5: (Color online). Qubit coherence (empty squares) in free- 
induction decay (FID), Hahn echo, and Carr-Purcell echo (CP) for a 
long "sinusoidal" spin chain with N = 1000 spins. The results from 
the 6-CCE (dotted lines) and 6-CE (solid lines) are compared with 
the exact solutions (empty squares). 



pairwise flip-flop between neighboring spins varies smoothly 
from one end to the other of the chain, so that a correlated 
cluster can grow to larger and larger size as time passes by. 
For a time greater than the bath spin flip-flop time r s f, the 
whole bath could become correlated. But the qubit decoher- 
ence would be completed within a time T much shorter than 
T s f , if the bath size is relatively large, or 

NT 2 r^ » 1, (33) 

due to the large number of small clusters contributing to the 
decoherence. Under the short-time condition, the CCE can be 
truncated with a rather small cut-off size M. The total contri- 
butions from clusters of various sizes to the decay of the qubit 
coherence as functions of the pulse delay time r [see Figs. [2] 
(b)-(d)] are shown in Fig. [4] For the free-induction decay in 
Fig. Ufa), the cluster contributions decrease rapidly with in- 
creasing the cluster size m at a time much shorter than t sI , so 
the 2-CCE (pair correlation approximation) already converges 
to the exact result. In the single-pulse Hahn echo [Fig.|4fb)], 
the higher-order correlations are more noticeable than in the 
free-induction decay, but the pair correlations still dominate 
for r <K r S f. For the two-pulse Carr-Purcell echo [Fig.[4fc)], 
as the decoherence due to the pair correlations is eliminated in 
the leading order of the spin-flip interactions, 19 the larger-size 
cluster correlations become important and a 6-CCE is required 
to reproduce the exact solution. 

In the relatively large bath, due to the large number of con- 
tributing clusters, the qubit decoherence completes when each 
individual cluster contribution is small, i.e., the small-term 
condition in Eq. d29l ) is satisfied. Thus the CCE coincides 
with the CE. This is verified in Fig. |5J where the exact solu- 
tion for qubit decoherence agrees with both the CCE and the 
CE truncated at the 6th order. 



i ' i ' r 




Pulse delay time 



FIG. 6: (Color online) Qubit coherence for a short "sinusoidal" spin 
chain of N = 100 spins in free-induction decay (FID), Hahn echo, 
and Carr-Purcell echo (CP). The results from (a) CCE and (b) CE 
truncated to the 2nd (dashed lines), the 4th (solid lines) and the 
6th (dotted lines) order are compared to the exact solution (empty 
squares). 



B. Small spin bath 

As the second example, we consider a short "sinusoidal" 
spin chain consisting of N - 100 spins. In this case the num- 
ber of contributing clusters is small and the decoherence pro- 
ceeds much slower. Fig. [6] shows that both the CCE and the 
CE converge to the exact result for a time much shorter than 
r s f. As the time approaches and goes beyond the bath spin 
flip-flop time r s f, the deviation from the exact solution is no- 
ticeable in both the CCE and the CE, indicating the emergence 
of correlations for clusters larger than the truncation cut-off 
size. But the CCE agrees with the exact result as long as it 
converges, while the CE converges to a different result (for 
T — t > 250 in free-induction decay, T — 2t > 300 in Hahn 
echo, and T — At > 800 in Carr-Purcell echo). The devia- 
tion of the converged CE from the exact result is due to the 
overlapping correction [see Eq. (1301 1 neglected in the CE: For 
T > Trf, the contribution of an individual cluster could be siz- 
able and the small-term condition in Eq. ( l29l > is violated, so 
the overlapping correction has to be taken into account. 

In the example discussed above, the CCE may not converge 
at a long time T > T s f. This is because the qubit-bath spin cou- 
pling assumes a smooth sinusoidal distribution and hence the 
energy cost of each neighboring pair is small and slow-varying 
as a function of the pair's position along the chain. The small 
and slow-varying energy cost of a pair-flip makes it possible 
for one pair-flip to affect its neighboring pairs, then the next 
neighbors, and so on and so forth. This way large-size cluster 
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FIG. 7: (Color online) Qubit coherence in free induction decay for a 
short "random" spin chain with N = 100 spins, (a) M-CCE and (b) 
M-CE are compared to the exact solution (empty squares). 
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FIG. 8: (Color online) The same as Fig.|7]except that the qubit-bath 
coupling constants jz„) are such that there are four neighboring spins 
in near-resonance. Notice that the 4th-order and the 6th-order trun- 
cations are not distinguishable by eyes in both CCE and CE. 



correlations may grow rapidly after the time goes beyond the 
pair-flip time r S f. In a relatively small bath of smooth qubit- 
bath coupling distribution, the CCE may fail to converge in 
the long time limit which is of interest in elongating the qubit 
coherence time by pulse control. 

The convergence problem of the CCE for small spin baths 
may be avoided if the qubit-bath coupling is random so that 
the bath correlation is localized and the size of correlated clus- 
ters is upper bounded. As the last example, we consider a 
short "random" spin chain consisting of N = 100 spins. In 
this case, the pair-flip energy cost is usually much larger than 
the spin-flip strength unless two neighboring spins are acci- 
dentally in near-resonance. Thus large-size cluster correla- 
tion can hardly grow significant even for a time well beyond 
the bath spin flip-flop timescale r s f. Figure [7] shows that the 
pair-correlation approximation already suffices for the qubit 
decoherence for an arbitrarily long time. Since higher-order 
correlations are small, the overlapping correction to the CE 
(which corresponds to at least three-spin cluster correlations) 
is unimportant, good agreement between the CE and the exact 
solution is also seen. Interestingly, in the small bath of ran- 
dom qubit-bath coupling, the qubit decoherence is not fully 
developed but coherent oscillations persist over a long time. 
These coherent qubit oscillations have indeed been observed 
in a system of a similar nature, namely, nitrogen-vacancy cen- 
ters in diamonds coupled to nuclear spins in the proximity. 40 
We would like to point out that the coherent oscillations can- 
not be reproduced by the LCE truncated as any finite interac- 
tion order. 

As discussed in Sec. Ill Dl the CCE and the CE differ in de- 



scribing the multi-spin correlations, which may be important 
for a small spin bath. Such difference can indeed be seen in 
Fig. |6]for a short "sinusoidal" spin chain. It is interesting to 
see the difference between the two theories in studying the co- 
herent multi-spin dynamics in a small bath. For this purpose, 
we consider a short (N = 100) "random" spin chain with acci- 
dental near-resonance between four neighboring spins. The 
pair-flip between the four spins in near-resonance may de- 
velop up to four-spin cluster correlations as the evolution time 
goes beyond the pair-flip timescale r s f. Indeed, Fig. Oa) 
shows coherent oscillations that are correctly reproduced by 
the CCE at the 4th or higher-order truncation. The correc- 
tion by the 5th and higher-order clusters is actually negligi- 
ble, verifying that the dominating cause of the coherent qubit 
oscillation is due to up to four-spin cluster dynamics in the 
bath. In contrast, the CE, though already converges at the 4th- 
order truncation, does not reproduce the exact solution [see 
Fig. Ob)], due to the neglect of the overlapping terms [see 
Eq. 63]. 



IV. CONCLUSION 

We have developed a CCE approach to solving the many- 
body dynamics of a generic interacting spin bath relevant to 
the center spin decoherence problem. In this approach, the 
bath propagator is factorized exactly into the product of clus- 
ter correlation terms, each of which accounts for the corre- 
lated flip-flops of a group of bath spins. In terms of the stan- 
dard LCE, a cluster correlation term corresponds to the infinite 
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summation of all the connected diagrams with all and only the 
spins in the cluster flip-flopped. For a finite-time evolution 
as in qubit decoherence, a convergent result can be obtained 
by truncating the expansion up to a certain cluster size. The 
CCE gives exact results whenever it converges. The lowest 
nontrivial order of the CCE corresponds to the previously de- 
veloped pair-correlation approximation. Compared to the CE 
method, the two theories yield similar results for large spin 
baths, but for small spin baths only the CCE accurately takes 
into account the multi-spin cluster correlations. As a simple 
method to sum over an infinite series of LCE diagrams, the 
CCE method can be readily applied to the case of interacting 
bosons or fermions. 
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APPENDIX A: LINKED CLUSTER EXPANSION FOR SPIN 
BATH DYNAMICS 

1. Thermal ensemble LCE 

The LCE for the propagator of an arbitrary spin bath in a 
thermal ensemble has been derived in Ref. 145}]. Here we 
summarize the main results including the Feynman diagram 
representation, as the basis of the extension to the LCE for an 
arbitrary noninteracting bath state. 

For a thermal ensemble of spin baths characterized by the 
noninteracting density matrix po = e^ H ° /Tie^ H " with Hq = 
Yj m u> m J z m , tne evolution in the contour time-ordered form is 



(Al) 



Here the contour time-dependent Hamiltonian H(t) switches 
alternatively between H <+ ^ and on the contour [see Fig. [2] 
(a)]. As an example, the pairwise bath Hamiltonian in Eq. (f2]i 
leads to 

H{t) = ^ + Yj Z nWt) + Y J D m , n {t)Jl(t)Jl{t) 



J] B m ^t)r m (f)j-{t), 



(A2) 



where Cl{t), Z„(t), D„, „(f), and B mA (t) depend on which con- 
tour time segment t is in, and the spin operators J z (t) = J z , 
J nit) = /* are time-independent but are written with explicit 
time-dependence to keep track of the time-ordering. 
Eq. ( IA1 b can be expanded into series as 



£ ens = Tr 



Po 



V fdti - - f dt„ %H{h) ■ ■ ■ H(t n 



(A3) 



where each term is the sum of ensemble-averaged ^-products 
of spin operators. Consider the ensemble average (F n ) = 



Tr [poF „] of an arbitrary product F„ = ^ c Jm\ ' " " ^m" (*«)> 
where the subscripts m q e { 1 , 2, • • ■ ,N} label bath spins and 
the superscripts a q e {z, +, — } label the spin operators. In or- 
der for {F„) not to vanish, J + and J operators must make 
pairs. An F n consisting of J z operators only is called a fully 
contracted product, whose ensemble average is trivially eval- 
uated. So next we consider the case with at least one J + oper- 
ator, 

For the moment, we assume that the product in F n is already 
time-ordered with t \ > ti > ■ ■ ■ > t„ and use the following pro- 
cedure to reduce (F„) to the sum of ensemble-averaged fully 
contracted products. First we move the spin raising operator 
JmStk) t0 tne r ig nt ' generating commutators between J^ t (tk) 
and the operators on its right, 



(F n ) = Tr 



Po 



2 JZ(h)--\j^k)jZ^)\--CStn) 



p=k+\ 



+ Tr [p J a m \ (*!)••■ (t k -i)J a Z\ fe+i) ' ' ' CS^JtM ■ 

Then in the second term, we use the cyclic invari- 
ance of the trace to move JL(tk) to the left (before 



po), and use /+(f)po 



e^'-poJ+it) to move JUt k ) 



across po back to its original position to get e^" m t (F„). 
During this process, we get additional commutators 
^ Tr [po Z k p l\ C\ (h) ■ ■ ■ [j+ k (t k ), C p (t p )\ ■ ■ ■ J%(t n )\ . Col- 
lecting all terms, we obtain 

(F„) = Tr po J] CSh) ■ ■ ■ [C p (t P )]' [J^(fk)]' ■ ■ ■ J%(tn) , 

- p(*k) 

(A4) 

where the contraction between a spin raising operator /*(f) 
and an arbitrary spin operator J"{t a ) is defined as 

[K(t tt )\' [J+(t)]' = Unita)]' 

with G m (t) = 0{t) [1 + f(to m )] + 8(-t)f(tD m ) being a Green's 
function and f(cj) = l/(e^ w - 1). Note that the contraction 
[J%{t a )]'[J+(t)] m ~ [J", Jm] (ta) is still an operator associated 
with a contour time t a and should be used in subsequent con- 
tractions. With [J Z ,J + ] = J + and [J~,J + ] = -27% the con- 
traction between /+(?) and J z m (t a ) [or J m (t a )] eliminates 
and converts J z m {t a ) [or J~(t a )] to [or -2J z m (t a )], reduc- 

ing the number of spin operators by one. The resulting product 
in Eq. (1A41 i is still in a time-ordered sequence, so the time- 
ordering operator can be recovered so that 

Tr [po&JZ) ft)-- - J«£ (tk-iVU^C^ fe+i) ' ' ' J %(tn)] 



= Tr 



po 2 (h) ■ ■ ■ [/,;»]• [jZ(t P )]' ■ ■ ■ j%(t n ) . 

(A5) 



In the above equation the time ordering t\ > ti > ■ ■ ■ > t n is no 
longer assumed, since all the spin operators commute in the 
,% -product. 
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FIG. 9: (a) Diagrammatic representation of the single-spin term 
Z„(t)J~(t), the diagonal interaction D m „(t)J^(t)J~(t), and the off- 
diagonal interaction B m „(t)J^(t)J~(t). (b) Second-order fully con- 
tracted diagrams. 



The contraction procedure in Eq. (IA5t can be repeated 
whenever there is still a spin raising operator left. Thus Wick's 
theorem follows: Under the average over a noninteracting 
thermal ensemble, a ^-product of spin operators can be re- 
placed by the sum of all possible fully contracted products 
which contains only J z operators. 

According to Wick's theorem, each term in Eq. (IA31 > gen- 
erates a series of fully contracted products. The fully con- 
tracted products can be visualized by Feynman diagrams with 
the following definition of constituent elements and construc- 
tion rules. 

1. A spin operator / + , J z is represented by a vertex 
as a filled circle, an empty circle, or an empty square, 
respectively; 

2. Each diagonal (off-diagonal) interaction term contain- 
ing n spin operators in the Hamiltonian is represented 
by a dashed (wavy) interaction line connecting n ver- 
tices; 

3. Each contraction [./^(f)]" [J%(ta)]' is represented by a 
solid arrow starting from the vertex J^if) and ending at 
the vertex J„{t a ). At the end of the propagating arrow 
the commutator [Jf n , /+] is to be taken; 

4. Each J + (f) vertex denoted by a filled circle is connected 
with one outgoing propagating arrow, each J z (t) vertex 
denoted by an empty square is either free-standing or 
connected to one incoming arrow [converting J z (t) to 
J + (t)] and one outgoing arrow [from the resulting oper- 
ator J + (t)], an d eacn J~(t) vertex denoted by an empty 
circle is connected to one incoming arrow [converting 
J~(f) to -2J z {t)~\ or two incoming arrows [the first ar- 
row converting J~(t) to -2J z (t), and the second arrow 
converting -2J z (f) to -2J + (t)] and one outgoing arrow 
[from the resulting J + (t)]. 

Note that each fully contracted product (and hence each dia- 
gram) is an operator consisting of J z spin operators only. 

Taking the Hamiltonian in Eq. (IA21 > for example, the 
diagonal single-spin term Z n (t)J z (t), the diagonal inter- 
action D mt „(t)J^(t)J^{t), and the off-diagonal interaction 
B mA {t)J^{f)J~{f) are visualized in Fig. |9ja). The first-order 
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FIG. 10: (Color online) (a) A third-order connected diagram consist- 
ing of one diagonal and two off-diagonal interaction terms, (b) A 
4th-order connected diagram consisting of four off-diagonal interac- 
tions. The contraction processes contained in each diagram are given 
below. 



expansion (-/) dt H(t) in Eq. jA3i gives two fully con- 
tracted products [the first two diagrams in Fig. |9ja)]. The 
second-order expansion of Eq. ( IA3I ) gives four fully con- 
tracted products shown in Fig.|9|b). The first three diagrams 

(i) = { ^z m (t x )ZM)J z m {h)W 2 ), 

{-if 

(ii) = 2 x K -^Z m (h)D n j(t2)J z m (h)J z „(t2)Jf(t2), 



(-0 



(iii) = ^A„,,^i)ZVfe)4('i)^('i)^te)^te), 

come from diagonal terms and involve no contractions. Here, 
as a convention, we have suppressed the sum over spin indices 
and the contour time integrals. The last diagram 

(-0 2 

(iv) =—rj-B mtn {t\)B m<n (t2) 

iWV [W]" [Km' [J-(t 2 )Y (A6) 

comes from two off-diagonal interaction terms and involves 
two contractions [J*(h)f [J~(t2)Y and [Jnih)]' Kfe)]*, a s 
indicated by the two solid arrows in Fig.|9](b) (iv). 

Similarly, higher-order diagrams can be constructed 
by using the above Feynman rules. Figure \W\ gives two 
examples. The 3rd-order diagram in Fig. [Tola) con- 
sists of two off-diagonal interactions B„ hn {ti)J^{t\)J m (_t\) 
and B„ un (t2)J^(t2)J^(t2) a nd one diagonal interaction 
Aafe)^n( f 3)-^( f 3)- K contains three contractions, one on 
spin J,„ and two on spin J„, corresponding to the three solid 
arrows. The 4th-order one in Fig. [TOl'b) consists of four 
off-diagonal interactions and contains five contractions, three 
on spin J„ and one on each of the other two spins. 

To illustrate the evaluation of the diagrams, we consider 
again the secular pair-interaction Hamiltonian in Eq. iA2\ as 
the example. The rules for constructing the analytical formula 
for a diagram are: 

1. A contour-time dependent constant is associated with 
each interaction line, namely, Z n (f) for an open-ended 
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FIG. 11: Topologically inequivalent connected diagrams up to tl 
4th order for the Hamiltonian in Eq. {2^. 



dashed line representing the spin splitting, D„ h „(t) for 
a dashed line connecting two diagonal spin operators 
representing the diagonal interaction, and B m „(t) for a 
wavy line representing the off-diagonal interaction; 

2. Each solid arrow from J„(t) to J" t (t a ) gives the Green's 
function G m (t a —f), each freestanding vertex gives J~, 
each J~ vertex connected to one incoming arrow gives 
(-2)J~, and each 7,7 vertex connected with two incom- 
ing arrows and one outgoing arrow gives (-2); 

3. A global factor {—i) k jk\ is associated with a diagram 
containing k interaction lines; 

4. The spin indices are summed over and the contour times 
are integrated over. 

For example, the last diagram in Fig.|9jb) gives 

(-0 2 



-£„,,„ (fi) B m ,„ (f 2 ) G,„ (t 2 -h)G n (h - ti) (~2J Z J (-2J z n ) , 



2! 



which can also be evaluated directly from Eq. (lA6t by carry- 
ing out the two contractions. 

Summation of all the fully contracted diagrams leads to the 
LCE of the ensemble-averaged evolution 
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FIG. 12: Mapping of spin-1 operators /* to pseudo-spin- 1/2 opera- 
tors. The so-defined mapping restricts the state evolution of the two 
pseudo-spin- 1/2's within the physical Hilbert space of the spin-1. 



2. Single-sample LCE for spin- 1/2 baths 

For a single noninteracting bath state \£T) (in contrast to the 
thermal ensemble), the LCE for spin- 1/2 baths has been out- 
lined in Ref . i23ll . Here we reproduce the main results using 
the ensemble LCE depicted in the previous subsection. 

The key is that for a spin- 1/2 bath, any noninteracting state 
\SF) = ®n \ jn) can be taken as the ground state of a correspond- 
ing noninteracting Hamiltonian, 



H 



(A8) 



(A7) 



where io„ < (or > 0) for \j n ) = |f) (or \\)). So the nonin- 
teracting single-sample average C7"|<9|J7") becomes the zero- 
temperature limit (j6 — > +oo) of the corresponding noninter- 
acting ensemble average Tr(pjO) with the density matrix 
pj = e~P Hj /Tr(e^ Hj ). In particular, the single-sample ex- 
pectation value of the bath propagator 

£=(jW e- i £ Hm \j)= lim TA Ps ^e- i ^ H(,)dt }. 

\ I I / (S— >+m L J 

Thus the single-sample LCE is obtained by simply setting the 
Green's function G„(t) - 9{t)5 - 6(—f)Sj n j and replacing 
the ensemble average Tr(po • ■ ■ ) with {JJ |- ■ ■ | J'), i.e., 

X = I J I %e-'' X H «<"| j\ = (j |«*| j) . (A9) 

Note that the connected diagrams in n contain only J z opera- 
tors, which commute with each other. Thus the single-sample 
average can be performed for each diagram to convert it into 
a c-number, so that 

i: = exp(CT|7r|J-» = exp(7r). 
3. Single-sample LCE for higher spins 



where n is the sum of all the connected diagrams, such as the 
first two diagrams in Fig. [2 a) and the last diagram in Fig.[9fb), 
but does not include the disconnected ones such as the first 
three diagrams in Fig.[9jb). As an example, for the Hamilto- 
nian in Eq. dA21 >. all the topologically inequivalent connected 
diagrams up to the 4th order are shown in Fig.Q~T] We see that 
the number of diagrams increases significantly with increas- 
ing perturbation order. 



For a higher-spin bath, a noninteracting single sample state 
|J7) in general is not the eigen state of a noninteracting Hamil- 
tonian as in Eq. ( 1A8I >. Thus the single-sample LCE for higher- 
spin baths may not be derived from the ensemble LCE di- 
rectly. Here we provide a solution by mapping a higher spin 
to a composite of pseudo-spin-l/2's. The mapping is such 
that the physical states form an invariant subspace in the much 
larger pseudo-spin Hilbert space. 
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FIG. 13: Diagram representation of (a) Z„J~, (b)D mn /^,/J, and (c) 
B m „J^J^. (d) Topologically inequivalent second-order connected di- 
agrams for a spin-1 bath. 



Without loss of generality, we consider a spin-1 J„. The 
mapping from the spin-1 states to the states of two pseudo- 
spin-l/2's S„j (A = 1,2) is 

|1,-1>„ -» U)„,il4>„,2, 

|i,o>„^IT) n ,ili>„,2, 
|i,+i>„ ^ |T>„,i IT>„, 2 , 

as schematically shown in Fig. [12] The two pseudo-spin- 1/2's 
have more basis states than the spin- 1 . To restrict the evolution 
of the spin state within the physical Hilbert space, we map the 
spin-1 operators to pseudo-spin- 1/2 operators as (see Fig.[T2"b 

£-» V2S: j (l/2-S-J + V2(l/2 + S; u )s„+ 2 , 
Jn - V2S- 1 (l/2-^ n2 )+ V2(l/2 + S^)s- 2 , 

J n ~~* S n,l + 5 n,2' 



in which the flip-flop of one pseudo-spin- 1/2 is conditioned on 
the state of the other pseudo-spin-1/2. Thus each spin-1 oper- 
ator is mapped to an interaction term of the two pseudo-spin- 
1/2's. An interaction term between two spin-l's would con- 
tain up to four pseudo-spin-1/2 operators, making the Hamil- 
tonian rather complicated. Nonetheless, the LCE can be read- 
ily applied to the spin-1 bath mapped to a pseudo-spin-1/2 
one. For a spin-1 Hamiltonian as in Eq. ( IA2I ). the different 
terms Z n (t)J z n (t),D n%ll (t)Jl n (t)J z n (t), and B„Jt)J*(t)J-(f) con- 
verted to pseudo-spin-1/2 operators, are represented in turn 
by the three diagrams in Fig. [Ol a). There each vertex (denot- 
ing a pseudo-spin-1/2 operator according to the same Feyn- 
mann rule depicted in Appendix I A II ) is associated with a 
physical spin index n e {1,2, •• • ,N] and a pseudo-spin in- 
dex A e {1,2}. The mapping leads to three different types of 
off-diagonal interactions (including simultaneous interaction 
involving up to four pseudo-spin-l/2's). As a result, the num- 
ber of connected diagrams increases dramatically, e.g., we 
have six topologically inequivalent second-order connected 
diagrams for spin-1 baths, as shown in Fig. fOl b), while for 
a spin- 1/2 bath we have only one such diagram (the third dia- 
gram in Fig.fTTT). 
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